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The stability of Bose gases near resonance has been a puzzling problem in recent years. In this 
Letter, we demonstrate that in addition to generating thermal pressure, thermal atoms enhance the 
repulsiveness of the scale-dependent interactions between condensed atoms due to renormalization 
effect and further stabilize the Bose gases. Consequently, we find that, as a precursor of instabil¬ 
ity, the compressibility develops an anomalous structure as a function of scattering length and is 
drastically reduced compared with the mean-field value. Furthermore, the density profile of a Bose 
gas in a harmonic trap is found to develop a flat top near the center. This is due to the anomalous 
behavior of compressibility and can be a potential smoking gun for probing such an effect. 

PACS numbers: 67.85.Jk, 67.10.Ba 


Recently, interest in interacting quantum gases has 
been revived due to the application of Feshbach reso¬ 
nance [H-Q. Feshbach resonance provides the tunability 
of the scattering length a, which uniquely characterizes 
the low-energy inter-particle interaction, from zero to in¬ 
finity and also from negative to positive. It offers an easy 
experimental access to resonant quantum gases. Since 
the conventional perturbation theories are no longer valid 
for these strongly interacting systems, we are confronted 
with the theoretical challenge of unravelling the puzzle 
of large-scattering-length physics. Among the various 
applications of Feshbach resonance, what has attracted 
particular interest is the atomic Bose gas at large posi¬ 
tive scattering lengths, known as the resonant Bose gas 
on the upper branch of a Feshbach resonance mi. 


For a dilute Bose gas where the scattering length 
is small and positive, the dominant contribution to 
its chemical potential is described by the Hartree-Fock 
mean-field value, (j-hf = 47ran, with n the number den¬ 
sity of the particles. The leading order correction to the 
mean-field result is proportional to a dimensionless pa- 


na^ <C 1. Near resonance where na^ ^ 1 or ^ 1, the 
dilute-gas theory [H-il is not applicable and several 
nonperturbative approaches have been taken to under¬ 
stand resonant many-body physics. They all predict the 
fermionization of bosons near resonance, which denotes 
that the chemical potential of the Bose gas reaches nearly 
the Fermi energy of a Fermi gas with the same density 
*281-^. In other words, the chemical potential of a Bose 


gas is predicted to increase linearly with a in the mean- 
field regime when na^ <C 1 and saturate at an energy 
scale defined by n as resonance is approached. These 
theory predictions so far are consistent with experimen¬ 
tal studies on chemical potentials Q. 


One of the major challenges of creating and probing 
unitary Bose gases experimentally comes from the short 
lifetime of the gases. Overcoming this difficulty has been 
one of the main focuses of a few recent experimental ef¬ 


forts [111, llJ, |18| . Especially, a very interesting attempt 
was made a year ago by quenching the scattering length 
and stabilizing the gas for relatively long time via a dy¬ 
namic approach 1^ . It had been generally believed that 
the short lifetime can be solely attributed to the few-body 
losses, which tend to increase rapidly when a resonance is 
approached. On the other hand, recent theoretical works 
suggest a more dramatic onset of many-body instability 
beyond a critical scattering length at zero temperature 


M l32| . This quantum critical point is mainly dictated 


by the scale dependence of interactions or the running 
of the coupling constant. Beyond the critical point, the 
coupling constant changes its sign at a relevant many- 
body energy scale, resulting in an effectively attractive 
interaction between condensed atoms and hence an in¬ 
stability. The competition between few-body losses and 
the onset of many-body instability remains to be further 
studied in future experiments. 


Given the limited current understanding of this issue, 
a few questions that are both fundamentally and practi¬ 
cally interesting arise. 1) What is the stability domain 
for Bose gases in the T — a plane? 2) What are the 
mechanisms for stabilization at finite temperatures? 3) 
What anomalous properties does a Bose gas have near 
the stability domain? 4) What can be used as a smok¬ 
ing gun for the stability boundary in experiments? In 
this Letter, we present an investigation of this topic and 
address the four questions raised above. Although there 
have been several previous theoretical studies on near¬ 
resonance Bose gases, the y w ere exclusively focused on 
the zero temperature case 28-3^. In contrast, the finite 
but low temperature case, which is more closely related 
to experiments, demands further investigation. 


Fig. d] shows the phase diagram for resonant Bose 
gases. The main plot displays the phase diagram near 
the quantum critical point (the red dot on the horizontal 
axis). The phase boundary is defined as the line where 
the compressibility changes its sign. It separates the 
regions with positive and negative compressibility. As 
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FIG. 1. (color online). The phase diagram of a scattering 
atomic Bose gas near the T = 0 qnantnm critical point. 
Temperature is shown in unit of the “Fermi temperature” 
Tf defined for a Bose gas with density n, which is equal to 
« 7.6n^^^. The blue solid line is where the 
compressibility k changes its sign. It separates the stable and 
unstable regions where the compressibility is positive and neg¬ 
ative respectively. The red dot on the horizontal axis repre¬ 
sents the quantum critical point. The inset shows the overall 
phase diagram in the same parameter space with both pos¬ 
itive and negative scattering lengths, where the red crosses 
represent the data from our calculation. The red dashed line 
connecting the data is drawn as a guide for the eyes in the 
vicinity of resonance. The dot-dashed box marks the corre¬ 
sponding region in the main plot. For reference, we plot in 
blue and black dotted lines the BEG transition temperature 
for an ideal Bose gas Tbec and Tp respectively. The yellow 
dotted line shows the position of resonance. 


shown in Fig. [1] the critical scattering length increases 
as the temperature rises from zero, indicating that small 
finite temperatures stabilize the gas. In our calculations, 
the compressibility k is calculated via, 


1 dn 
K dn 


( 1 ) 


from the equation of state p.(n). /x is obtained from 
= dF{nQ, ^)/dnQ, where Uq is the number density of 
condensed atoms and F{nQ,^) is the free-energy density 
F{no,fi) itself is calculated at a prefixed chemi¬ 
cal potential, thus is determined self-consistently (refer 
to Eq. ([U and the context there). In the low-temperature 
limit, we have identified two dominant processes that 
yield the leading order contributions to F. We denote 
them as F Ri Fi -b F 2 . Fi has the following form, 


J ( 2 ) 


where tk = fc^/2 is the energy of a free particle, rj = E—/i, 
E is the self energy, and (3 = 1/T. For clarity, we set the 


reduced Planck constant h, the Boltzmann constant 
and the atomic mass m to be unity throughout the paper. 
Fi counts the contributions of scattering thermal atoms 
whose energy receives a modification of 77 due to many- 
body interactions. As expected, the thermal pressure 
provided by Fi enhances the thermodynamical stability 
of the Bose gas. F 2 on the other hand can be attributed 
to the interaction energy of the condensed atoms and can 
be written as 

F 2 = ]^92{ri)nl, (3) 

where 92 ( 11 ) corresponds to the renormalized two-body 
running coupling constant at the energy scale defined by 
rj [ 3 ^. At finite temperatures, 92 ( 9 ) has the form 


92 iv) = 


Aira 


47r 


d^q 2nB(q^/2 + r)) 

(27r)3 q^ + 2r] ’ ^ ’ 


where nsix) = (e^“ — is the bosonic distribution 
function. The last term of Eq. (jd]) represents the bosonic- 
enhancement effect due to thermally excited atoms. This 
makes the two-body interaction more repulsive and in 
turn stabilizes the gas. In the limit of low temperature, 
we also perform an analytical expansion of Eqs. m and 
©■ Our analysis shows that the thermal pressure is the 
driving force that stabilizes the gas only when T is small. 
At r Ri O.ITTf, the contribution of the running coupling 
constant becomes comparable to that of thermal pressure 
and hence more dominant. 

When the scattering length is tuned to be negative, 
the interaction between atoms is effectively attractive. 
At low temperatures, the gas is thermodynamically un¬ 
stable. It has been shown that for a small negative scat¬ 
tering length, Bose gases become stable at a tempera¬ 
ture slightly above the transition temperature of Bose- 
Einstein condensate (BEG) for an ideal gas, Tbec |36- 
. This result can also be obtained using our approach 
outlined here. When T > TbeCi n is related to F by 
n = —dF/dfj,. By further supplementing an equation 
for the self energy E = Siran, we can identify the phase 
boundary for small negative scattering lengths. Now we 
present a more thorough phase diagram for Bose gases, 
shown in the inset of Fig. [TJ The red crosses mark the 
instability points where k changes its sign for both small 
negative scattering lengths and also around the quantum 
critical point. Near resonance, we extrapolate between 
these two cases by smoothly connecting the two sides. 
The resultant curve (red dashed line) near a = 00 is qual¬ 
itatively consistent with the picture in Ref. derived 
via a high-temperature expansion. 

Now we turn to questions 3) and 4) asked in the in¬ 
troduction regarding the experimental signature of the 
stability boundary. A previous study on two dimensional 
Bose gases which applied a theoretical approach simi¬ 
lar to the one outlined here implies an anomalous behav¬ 
ior of the compressibility. That is the inverse of com¬ 
pressibility as a function of scattering length develops 
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FIG. 2. (color online). The change of sign of compress¬ 
ibility as functions of (a) scattering length when tempera¬ 
ture T — 0.2Tf and (b) temperature when scattering length 
a = 0.1739n“^^®. (c) Compressibility as a function of density 
for fixed scatteirng length at temperature T — (blue 

solid line). The compressibility drops rapidly near the critical 
na^, compared with the result from Lee-Huang-Yang (LHY) 
theory [T^, (red dashed line). 



FIG. 3. (color online). Density profile for N = 9.06 x 
atoms in a harmonic trap with frequency tu = 3 x 
The blue solid line shows the result of our calculation within 
local-density approximation at temperature T = 

The density in the trap center is no? = 5.25 x 10“^, slightly 
below the critical value. Note that the rapid drop of com¬ 
pressibility near na® = 5.25 x 10“^ (see Fig. [2](c)) leads to 
the flat top of the density profile. We also plot the density 
profile from LHY theory in red dashed line for comparison. 


a maximum as a precursor of the instability. This is 
consistent with a recent experiment in Chicago on two 
dimensional Bose gases [l^ where a maximum in the in¬ 
verse of compressibility was measured. Given the success 
in two dimension, here we further show in Fig. the 
data for the anomalous compressibility along the phase 
boundary. Figs. [D (a) and (b) show the sign change of 
the compressibility. Before reaching the instability point, 
the system has a positive compressibility and is thermo¬ 
dynamically stable. As a precursor of the instability, the 
compressibility decreases rapidly as shown in Fig.[^(c). 

Furthermore, since the gas with a negative compress¬ 
ibility is difficult to control, it is appealing to ask what 


is the potential smoking gun for the stability boundary 
in a three dimensional Bose gas before it becomes unsta¬ 
ble. For trapped Bose gases, this anomalous compress¬ 
ibility leads to a distinct feature in the density profile. 
Within the local-density approximation the trap renders 
a position-dependent chemical potential, which reaches 
its maximum value in the center of the trap and decreases 
towards the edge [i^. Consequently, the density profile 
of the gas, n(r), is a direct reflection of the compressibil¬ 
ity and can provide experimental evidence of the quan¬ 
tum critical behavior discussed above. In Fig. [31 we plot 
the density profile for a Bose gas in a harmonic trap. The 
peak density, which is located in the center of the trap, 
is set slightly below the critical value na^ = 5.25 x 10“^. 
N = 9.06 x 10"^ is the critical particle number for this set 
of to (the trap frequency), T, and a, beyond which the 
gas near the trap center becomes unstable. In general, 
we find that N^r ~ f {Ta^), where Ncr is the 

critical particle number and / is a dimensionless func¬ 
tion and /(O) = 0.0024. For comparison, we also plot 
in the red dashed line the density profile obtained using 
the Lee-Huang-Yang (LHY) theory [3, [13 for the same 
w, a, and N at zero temperature. It can be seen that 
our calculation predicts a much flatter top for the den¬ 
sity profile compared with the LHY result. This is due 
to the rapid drop of k near the critical point (i.e., the 
gas becomes difficult to compress before reaching the in¬ 
stability), in contrast to the mild decrease in the LHY 
theory, as shown in Fig. (c). There have been exper¬ 
imental evidences of a flat top of the density profile for 
a Bose gas in a harmonic trap near resonance, but they 
have long been attributed to the strong three-body loss in 
the center of the trap. However, our analysis shows that 
the worrying three-body loss is still relatively small near 
the instability point, which we will discuss in more de¬ 
tails below. As a result, the observed flat top is a strong 
evidence of the anomalous compressibility, a precursor of 
the onset of many-body instability. 

Bose gases of scattering atoms suffer three-body- 
recombination processes due to the presence of bound 
states deeper in the energy spectrum. When three atoms 
scatter, two atoms can form a bound state while the third 
atom gains appreciable kinetic energy and may escape 
from the trap; this is known as three-body loss. The gas 
therefore has a finite lifetime. The lifetime due to losses is 
determined by the rate of this three-body-recombination 
event, known as the three-body loss rate L 3 . There have 
been extensive studies on L 3 [ill, [3 41-43|. It was shown 
that, L 3 increases as (apart from the peaks and dips 
due to interference effects) and saturates to ^ when 
a becomes comparable to the thermal wave length. This 
smooth variation of L 3 with a is fundamentally different 
from the sharp sign change of k at the onset of the many- 
body instability discussed above. In order to further com¬ 
pare the few- and many-body effects, we carry out an 
estimation of L 3 near the quantum critical point. In the 
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low-temperature limit, L 3 was shown to be 203.7a^ up to 
a periodic function of In(aA), oscillating between 0 and 1 , 
where A is an ultraviolet momentum scale depending on 
the short-distance details of the interatomic interaction 
M\. Near the critical point where ~ 0.174, the en- 
ergyscale set by the three-body loss rate can be estimated 
as L^n'^ Ki 0.024:6Tp, while the chemical potential near 
this point is calculated to be /r ~ 0.735 Tf, i.e., the onset 
of many-body instability sets in before the three-body 
loss becomes significant. This is in contrast to the com¬ 
mon belief that few-body loss is the main limitation to 
the experimental accessibility of highly degenerate Bose 
gases. On the other hand, we can also infer that, com¬ 
pared with the many-body timescale, the lifetime due to 
three-body loss is long near the critical point. This vali¬ 
dates our treatment of the system as a thermodynamical 
system and justifies our thermodynamical analysis above. 
Besides, the presence of the ultraviolet parameter A in L 3 
reflects the nonuniversality of Bose gases near resonance, 
which is related to the Efimov physics [i^. This effect 
can be incorporated in the effective-field-theory approach 
quite straightforwardly at zero temperature by including 
the three-body interaction energy with an ultraviolet cut¬ 
off [3l|. By varying the ultraviolet cutoff from 
to our calculation shows an oscillation of the 

position of the critical point, indicating nonuniversality, 
but the relative magnitude of the oscillation is less than 
2%. This suggests that the three-body effect is negligi¬ 
ble even near the critical point, which is consistent with 
previous studies. 

At last, we briefly introduce our technique. We adopt 
an effective-field-theory app roach combined with self- 
consistent equations [3l|, l32| . In this approach, the free- 
energy density F{no, fj,) (which reduces to the energy den¬ 
sity at zero temperature) is calculated as a function of 
condensed-particle density hq and a preassumed chemi¬ 
cal potential fi. At equilibrium, fj, is further related to 
F{no,n) by 


dFino,^) dF{no,n) 


where the second equation is the number equation. The 
chemical potential, condensation fraction, and thermo¬ 
dynamical quantities like compressibility can be obtained 
by solving Eq. ([5]) self-consistently. 

The evaluation of F can be carried out diagrammati- 
cally [Hillmill, 


F{no,fi) = 


00 (0)/ \ 

9 m W 

^ M\ 

M=0 


n, 


M 

0 J 


( 6 ) 


where effective M-body interaction. How¬ 

ever, this traditional perturbation expansion converges 
quite slowly. In our approach, we further introduce the ir¬ 
reducible A/-body interactions which depend on uq them¬ 


selves [3l|. In this case, F can be rewritten as. 


00 ^(IR) 

F{nQ,g) = ^ 

M=0 


9m {9,no) m 

I Li 


M\ 


0 > 


(7) 


/TD \ 

where is the irreducible M-body interaction. At 
zero temperature, this summation starts with M = 2, 
and 52 ^^^ has been shown to produce more than 99% of 
the LHY correction in the dilute limit. Near resonance, 
this expansion has also been proven to converge rapidly 
IH - According to our calculations at zero temperature, 
the contribution of g^^^^ to F near the critical point is 
less than 5% of that of ( 72 ^^^. A rigorous solution at 4 — e 
dimensions further demonstrated that the effect of 
M-body (M > 3) interactions is indeed suppressed by an 
extra power of e, which, after extrapolation, is consistent 
with the numerical smallness we observed here. In our 

present analysis, we neglect contributions with M > 3. 

fIR.) 

On the other hand, there are contributions from gQ at 
finite temperatures, which correspond to the free energy 
of thermal atoms. Among this type of diagrams, those 
involving a single thermal atom give the leading-order 
contribution while others are suppressed by additional 
factors. For example, the Nozieres-Schmitt-Rink 
type diagrams considered in Ref. [H| are suppressed in 
the low-temperature limit. 

In summary, we have studied the compressibility and 
stability of the resonant atomic Bose gas. Near the insta¬ 
bility line, besides generating a thermal pressure, thermal 
excitations further stabilize the gas by enhancing the re¬ 
pulsiveness of the interatomic interactions. Unlike the 
smooth varying of lifetime due to few-body loss, this on¬ 
set of instability due to the many-body effects is charac¬ 
terized by a sharp sign change of the compressibility. As 
a precursor of the instability, the compressibility drops 
quickly, which induces a flat top in the density profile in 
the presence of a harmonic trap. This might be a po¬ 
tential smoking gun for the onset of many-body instabil¬ 
ity. It is worth emphasizing that the instability studied 
here occurs when the few-body losses are still insignifi¬ 
cant; this makes it plausible to experimentally detect this 
many-body instability. Our study provides useful insight 
on the thermodynamics of scattering atomic Bose gases 
for a wide range of scattering lengths and can further 
shed light on future research on the dynamics of Bose 
gases. 
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